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Abstract. Sensitive measurements of the linearly polarized spectra of stars can be 
used to deduce geometric properties of their otherwise unresolved circumstellar envi- 
ronments. This paper describes some of the evidence for optical pumping and absorp- 
tive linear polarization and explores some interesting applications of linear spectropo- 
larimetry for obtaining spatial information from imbedded stars. 



1. Introduction 

Linear spectropolarimetry of stars obtained with polarization sensitivity of 0. 1 % or bet- 
ter of t en shows linear polarization variabi l ity acro ss absorption lines (cf. IVink et al7 " 



l2002l l2005bl : lHarrington and Kuhn I [20071 l2009al lbh. This effect is observed over a 
broad range of stel lar conditions and in many lines in a stellar spectrum. We've shown 
( Kuhn et al. 1120 07) that the anisotropy of the stellar radiation field can lead to anistropy 
in the ground-state atomic energy levels responsible for the observed stellar absorption 
lines. This qua ntum mechanical effect is called optical pumping (OP) in the laboratory 



dHapper 1119721) . and has been know n for some time to yield observable linear polar - 
ization effects in the solar spectrum dTrujillo Bueno and Landi DeglTnnocenti II 19971) . 
Using detailed radiative models of Hydrogen absorption above the surface of a star 



we've shown (Ku hn et al. 1 12007) that H a is often quite sensitive to optical pumping. 
Our earlier models didn't include the effects of hyperfine coupling and considered only 
the low-lying atomic levels in Hydrogen. 

Interpreting H a polarization is complicated by the degenerate hyperfine transitions 
that are not resolved in typical stellar spectropolarimetry measurements. In this paper 
we'll outline possibilities for simplifying the interpretation of the observed polarization 
across stellar absorption lines. For example, from certain non-hydrogenic lines we can 
show that optical pumping, not scattering polarization, is the mechanism responsible 
for producing spectropolarimetric variability. 

In stellar systems where the lines in the B aimer or Paschen sequence originate 
from the same circumstellar region, we find further support for the OP model by show- 
ing that the linear polarization ratios for the strongest hydrogen lines can be computed 
from simple models. These polarization ratios may be used to constrain the local cir- 
cumstellar environment. We examine several stellar systems where optical pumping in 
the B aimer lines appears to originate from the same region and we speculate on how 
detailed comparison of optical pumping in different lines can be developed as a useful 
tool. 
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In stellar systems with a rotating disk we illustrate how the absorptive linear po- 
larization due to optical pumping can be used to determine some basic geometric prop- 
erties. Thus, the orientation of a disk in the sky and its tilt with respect to the plane of 
the sky can be inferred from the absorptive linear polarization line profiles. Finally we 
show how linear spectropolarimetric data can be inverted to obtain spatially resolved 
"images" of the intervening gas near the imbedded star. 

2. Absorptive Polarization 




Figure 1. Optical pumping geometry for a cloud (blue ellipse) observed along 
a line-of-sight (red arrow) and radiatively excited from anisotropic stellar radiation 
(white arrow) 

Optical pumping originates from the lower level of atomic transition levels where 
the equilibrium density of the magnetic quantum number states depends on the mag- 
netic quantum number m. For example, if a lower-level state with J= 1 has a different 
electron density for the m = 0, compared to the m = -1, +1 levels then these atoms 
will yield a line opacity which depends on the incident light direction and linear po- 
larization states. The population anisotropy is determined by the angular anisotropy 
of the incident radiation. The optical pumping (OP) mechanism geometry we describe 
can be simplified to a "pumping beam (or light cone) represented by the mean direction 
of the radiation which induces anisotropy in the atomic sublevels, and a "probe" beam 
which is the direction of the optical radiation that is selectively absorbed by the pumped 
atomic gas in the line-of-sight to the observer. 

Figure Q] shows schematically how OP works in the circumstellar environment of 
a star. A classical analogy is indicated in this figure. We consider an absorptive "cloud" 
between the observer and a star indicated by the blue ellipse. The pumping radiation 
originates along the direction of the white arrow and effectively polarizes the dipole 
response of the atomic system to lie along this direction, indicated schematically here 
by the elongated blue ellipse. Light which is absorbed by the cloud originates from 
the photosphere of the star along the line-of-sight direction to the observer (the red ar- 
row). The absorbed radiation induces an electronic dipole response in the cloud which 
favors a direction along the elliptical long axis. Thus, for initially unpolarized radia- 
tion it is the probe beam electric field direction most closely aligned with the ellipse 
which is selectively absorbed and scattered. The perpendicular electric field orientation 
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is transmitted with lower opacity. The red ellipse projected in the side panel shows 
the enhanced polarization direction of the transmitted light as seen from the observers 
perspective. In general OP absorptive polarization yields a transmitted beam polarized 
perpendicular to the plane containing the pumping and probe beam directions, although 
a negative polarization can be induced in some atomic transitions as we show below. 



2.1. Database 

Our group has amassed a considerable database of precision linear spectropolarime- 
try measurements from a broad range (about 100) of "imbedd ed" stars. This includes 
HAeBe's, Po st-AGB, RV-Tau's, Wolf-Rayet, and others dHarrington and Kuhn ||2007 . 



I2008L 12009 al ibi). Part of this database has been obtained with our own HiVIS spectro- 
graph on the Haleakala AEOS 3.7m telescope but it also includes data obtained from the 
ESPaDOnS instrument on the CFHT. The data described below was obtained with the 
CFHT and comes primarily from archival data records. These data were obtained with 
spectral resolution R>50000 and with typical 2cr linear spectropolarimetric accuracy of 
0.1% in the line profiles. 



2.2. Non-scattering Lines 

Various mechanisms have been proposed to describe the linear spectropolarimetric vari- 

ability seen in stellar absorption lines (cf. Vin k et al. 120051 : [McLean Il979k l0udmaiier and Drew 



Il999l) . For example, sc attering polariza tion from free electrons or resonant line scat 
tering is often invoked dMcLean Ill979h . Unpolarized line scattering in combination 



with pola rized Thomson scattering can also yield polarization variability across spec- 
tral lines dMcLean II 19791) . Scattering polarization requires an atomic upper level which 
allows for a population anisotropy. This occurs if the net angular momentum (including 
nuclear spin coupling) yields a total angular momentum greater than 1/2. Calcium is 
relatively abundant with no net nuclear spin. The Call triplet transitions at 849.8, 854.2, 
and 866.2nm are an excellent multiplet to test scattering as a polarization mechanism 
since the lines have similar for mation properties and the longest wavele ngth term can- 
not be polarized by scattering. Ma nso Sainz and Trujillo Bueno I ((2003) studied these 
transitions in the solar spectrum in order to explain observed solar spectral polarization. 

The 866.2nm transition upper level is a J=l/2, 3 P\/2 state, while the other transi- 
tions originate from the P3/2 level. The lower level in each of the triplets is a J=3/2 or 
5/2 2 D level. Thus if line polarization is detected in the 866nm line it cannot be due to 
scattering polarization. The angular momentum of the lower levels in all of these tran- 
sitions allows for optical pumping. Thus, by detecting linear polarization in the reddest 
of the Call near-IR triplet lines, we have clear evidence of OP absorptive polarization 
which cannot be caused by resonant scattering effects. 



2.3. Call NIR Triplet results 

The Call triplet lines in many of our stars are in emission, but of those systems with 
these lines in absorption, about half show complex Stokes Q,U polarization structure. 
Invariably the linear polarization morphology across the red triplet component is similar 
to the two bluer triplet components. An example is illustrated here. U Mon is an 
RVTau 9 Id variable star. Fig. |2]shows the intensity spectra and continuum normalized 
polarized spectra. 
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Figure 2. Spectra and continuum normalized Stokes Q and U spectra for U Mon 
in the Call triplet lines 



The Stokes Q and U polarized spectra show what we now know to be typical 
spectral line polarization variation, localized in the absorptive part of (in this case) the 
Call triplet line cores. The linear polarization exhibits a typical amplitude of a few 0. l's 
of percent. The correspondence between the 866.2nm and the bluer triplet polarized 
spectra is clear evidence that it is an optical pumping mechanism which is responsible 
for the observed absorptive polarization signal. 



2.4. Optical Pumping and Hyperfine Levels in Hydrogen 

Hyperfine coupling of the atomic nuclear spin to the electronic levels implies that even 
J=0 ground-state levels may be sensitive to anisotropy in the radiation environment of 
a star. We've included this effect in hydrogen and built a radiation coupled model that 
includes 364 level s and about 5000 transitions up to n - 6 in hydrogen. Following 



(IKuhn et al. 1120071) we've solved for the m-dependent equilibrium population in each 
of these levels as a function of the background radiation field anisotropy (parametrized 
by the radiation "cone-angle") and the radiation temperature. We assume that atomic 
collisions and magnetic field quantum coherence effects can be neglected. For the suffi- 
ciently intense radiation field and low density stellar environments that we consider here 
these are generally good assumptions. Thus Fig [3] shows that a 10,000K black-body 
with a 73deg illumination angle produces higher populations in the largest m-states at 
larger L. Note that the n=l, L=0 ground-state has opposite m-dependent anisotropy, 
that is, the m=0 populations are greater than higher magnitude m states. The rela- 
tive population of higher n and L states is small but for more intense radiation (higher 
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temperature or larger illumination angle) the larger n and L states are relatively more 
heavily populated. As the cone angle (isotropy) increases the m-dependent population 
anisotropy decreases. 
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Figure 3. The relative population of hydrogen levels for a 10,000K black-body 
illuminating a collisionless gas with a 73deg cone angle. Hyperfine coupling is in- 
cluded and the curved lines show the relative population from -m to -Hm where m is 
the maximum magnetic total spin quantum number. The plot shows this population 
versus principal quantum number n and orbital angular momentum L. Thus for each 
n,L there are four levels with F=L+1, 2xL, and L-l. The thicker horizontal lines un- 
der each group of 4 total spin levels shows the relative population versus L for given 
n. The thick vertical lines are proportional to the relative population of all levels 
at given quantum number n. For example, here a 10000K black-body leaves most 
of the electrons in n=l states. Only the L=0,1 states have significant population at 
larger n. 



A strong m-dependent density dependence in radiation coupled levels in an atom 
can yield opposite population anisotropy in the upper and lower transition states. This 
is clear in Fig.[3]where the low lying L=0 states have higher m=0 populations, while the 
L= 1 states have lower m=0 populations relative to higher magnitude magnetic quantum 
numbers. In this figure a positive or negative anisotropy is indicated by states with 
line curvature up versus down. A consequence of this is that as the radiative coupling 
changes because of a stronger radiation field (pumping region is closer to the star), or 
hotter radiation, the OP polarization changes rapidly. Thus, changing anisotropy (as 
radiation illumination cone angle changes) affects a competition between the opposite 
anisotropy composite states to produce significant net OP line polarization changes. 

Polarization ratios of lines from the same term can be a very sensitive probe of 
the radiation strength and isotropy due to the opposite density anistropy that is possible 
in low-lying atomic states. For example, we have computed the ratio of the theoretical 
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Balmer line polarization in tLj/H (! ,. Figure |4] plots the ratio of the n=4 to n=2 and the 
n=3 to n=2 transition polarizations versus the radiation illumination angle in degrees 
for a 10000K blackbody. The dotted curve shows that while the line intensity ratio is 
nearly constant as the radiation isotropy and strength vary, the polarization ratio even 
changes sign in a complex way. Line polarization ratios are a sensitive probe of the 
near stellar environment and contain radiation information distinct from the intensity 
line ratios. 



o 



1 r 



r 




20 40 60 

Radiation angle [deg] 



50 



00 



Figure 4. The solid line shows the ratio of H^/Hq, polarization and the dotted line 
shows the corresponding line intensity ratio versus radiation illumination cone an- 
gle. Competition between opposite density anisotropy low-lying hydrogenic levels 
causes a complicated polarization ratio as radiation coupling varies with cone angle. 



Two examples illustrating a range of polarization and line ratios are SS Lep and 
MWC 361. SS Lep is a symbiotic binary star and MWC361 is a Herbig Be star with 
a directly detected disk (|Okamoto et al. Figure [5] shows H a and EL; for these 

two stars. The intensity ratio of the absorptive component of SS Lep and MWC 361 
are about 0.3 and 0.2, while the corresponding polarization ratios are 0.3 and 0.1. Our 
work to uniquely interpret such line and polarization ratios is only beginning. 



2.5. Geometric Properties from Absorptive Linear Spectropolarimetry 

The mean radiation pumping direction at any point in the circumstellar cloud defines 
a direction projected onto the plane of the sky that, in the simplest approximation (see 
Fig. [B, is perpendicular to the mean electric field polarization of the transmitted light 
through the cloud at that location. If each location in the cloud encodes a definite 
doppler shift velocity, for example if the cloud is actually a rotating disk projected 
against the stellar disk, then the shape of the polarized absoiptive line spectrum with 
wavelength gives information on the position angle of the disk in the sky and its incli- 
nation angle toward the observer. 

Figure[6]shows an example of the expected Q, U polarization for a disk oriented as 
in the figure and observations from e Aur that are consistent with this orientation. Sim- 
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Figure 5. Each panel shows the Q, U, and I spectra near H ff and forMWC 361 
and SS Lep. 



ilarly Figure |7] show the likely orientation of 51 Oph derived from H ff linear absorptive 
spectropolarimetry. 

For a star with sufficient symmetry in its disk a straightforward procedure to de- 
rive the disk orientation is to plot the polarization angle at each wavelength across the 
absorptive line profile. Thus, in the case of MWC 361, if we derive an angle from 
9 = l/2arctan(Q/U) and a polarization amplitude from p = ^Q 2 + U 2 then a plot 
of p(9) looks like Figure [8] The orientation of the disk is found by rotating the plot 
until symmetry about the x-axis is achieved. The data in Fig.[8]have been rotated about 
5 degrees to yield symmetry around the x-axis suggesting that the position angle of 
the MWC 361 disk is 5 d eg. This angle is consistent with the extended IR disk seen 
bv lOkamotoetal. I(l2009h but is at odds with the modeled rotation axis obtained from 
Zeeman Doppler Imaging (cf. lOkamoto et al."l l2009). 



2.6. Absorptive Linear Polarization Imaging 

With sufficient velocity resolution across the absorptive line profile the intensity and 
polarization amplitude and Q, U spectra may be inverted to yield the column density 
projected onto the star. Because the linear polarization is localized to the line-of-sight 
to the star, and analogous to the case of Doppler Zeeman Imaging for a stellar magnetic 
field, there is spatial information about the structure of the obscuring gas in front of the 
star in the Q, U spectra. In general the OP signal originates close to the pumping source 
- Kuhn et al. (2007) argued this is likely to be within 1 stellar radius of the photosphere. 
Either from the polarization amplitude or the absorbed intensity at each wavelength we 
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Figure 6. The observed Q and U spectra for e Aur and a modeled linear polariza- 
tion line profile for the indicated disk geometry. The small red lines in the disk are 
perpendicular to the transmitted light linear polarization direction at that projected 
disk location. 

derive a local cloud "density" or optical depth. At each wavelength the Q, U polariza- 
tion ratio defines an angle on the disk and the doppler velocity defines a surface (based 
on an assumed 3-space rotation profile or a wind velocity) the obscuring material must 
lie along. Alternatively the local polarization amplitude yields the separation between 
the local cloud position and the geometric center of the illuminating photosphere. We 
can demonstrate this using the H a spectropolarimetry for AB Aur in order to obtain an 
"image" of the intervening cloud as shown in Figure |9l Here a radial outflow wind ve- 
locity was assumed and the inversion has a 1 80 degree ambiguity because of the linear 
polarization 2-fold rotational symmetry. The light shades in this figure indicate where 
the largest absorption against the stellar disk occurs as obtained from a single epoch 
inversion using high resolution I, Q, and U H ff spectra of AB Auriga. 

3. Conclusions 

Absorptive linear polarization is a nearly ubiquitous phenomenon in the spectra of 
imbedded stars. It is often visible in polarized Stokes Q and U spectra when a sensitivity 
better than 1% is achieved. We've shown that it is not a scattering phenomenon and is 
well described by optical pumping due to ground-state anisotropy from the anisotropic 
stellar radiation environment. We've shown here that there is new information on the 
stellar environment contained in linear spectropolarimetry that is distinct from, for ex- 
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Figure 7. The observed Q and U spectra for the star 5 1 Oph and a correspondingly 
modeled linear polarization line profile for the indicated disk geometry. The small 
red lines in the disk/ring figure panel are perpendicular to the transmitted light linear 
polarization direction at that location. 



ample, line intensity ratios. With simple assumptions, the geometry of stellar disks 
can be inferred. With additional assumptions we've demonstrated how high resolution 
spectropolarimetric data can be inverted to yield "images" of the intervening gas from 
polarized absorptive line profiles. 
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Figure 8. The H ff and Hg Q, U amplitude and direction are plotted as a collection 
of points across the absorptive line profile. Rotating these points to a symmetry line 
around the x-axis shows that the disk position angle is about 5 degree in the sky. 
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Figure 9. From the high resolution I, Q, and U spectra of H a we derive the pro- 
jected optical depth of intervening hydrogen gas in the side panel of this figure. 
Lighter shades indicate higher depth. 



